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Abstract 

By generalizing the algebra of operators of the Asymmetric Simple Exclusion Process 
(ASEP), a multi-species ASEP in which particles can overtake each other,is defined on 
both open and closed one dimensional chains. On the ring the steady state and the 
correlation functions are obtained exactly. The relation to particle hopping models of 
traffic and the possibility of shock waves in open systems is discussed. The effect of the 
boundary condition on the steady state properties of the bulk is studied. 
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1 Introduction 



The aim of this letter is to introduce a natural multi-species generalization of the one- 
species asymmetric simple exclusion process (1-ASEP) [1-5]. We will show that using 
the Matrix Product Ansatz (MPA) for one dimensional stochastic systems, one can al- 
gebraicly define a natural p-ASEP ( p = number of species) such that all the results of 
1-ASEP will be reproduced in the special case p — 1. We will also show that this model 
incorporates new features which compared with 1-ASEP makes it a better candidate as a 
model of traffic flow. First we resume the basic facts about 1-ASEP. In this process which 
is defined on a one dimensional lattice with either periodic or open boundary conditions 
, particles hop stochastically to their right nearest neighbor site if this site is empty 
and stop if this site is occupied. This model is related to many interesting physical prob- 
lems including interfacial growth [6], directed polymers in random media [7,8], kinetics of 
biopolymerization [9], and traffic flow [10-12] or any queing problem. It is an example of 
a system far from equilibrium for which many exact results are known (for recent reviews 
see [5, 13]). Some of the basic results are as follows. 

On a ring of N sites where there are M particles interacting with each other, the steady 
state is such that all configurations have equal weights [14]. All the equal time correlation 
functions are also easy to calculate [14], that is 



Here the stochastic variable is if site k is empty and 1 otherwise. Therefore < r k > 
is the average density and < t^t; > is the density-density correlation. The unequal 
time correlation functions are not known in general and only certain quantities, like the 
diffusion constants have been calculated exactly [15,16]. 

On an open chain where particles arrive at the left with rate a and leave the right with 
rate (3 the steady state and the phase diagram [17-19] is known exactly. Particularly 
when a + f3 — 1, the steady state is defined by a Bernouli measure. The average density 
of particles is given by a and the steady current is J = cc(l — a). The measure being 
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< TiTj > = 
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N(N -1) 
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multiplicative, there is no correlation in this case. 

When a + (3 7^ 1, the steady state is again multiplicative [19 ], but this time matrices 
replace numbers. This is the content of the matrix product ansatz (MPA)[20]. In this 
case the probability of a configuration (ti,t 2 , ...,t n ) is given by 

P(n, T 2 , ..,7*) = < W\D{t 1 )D{t 2 )...D{t n )\V > (1) 

where the matrices D := D(l) and E := D(0) and the vectors < W\ and \V > satisfy 
the following relations: 

DE = D + E D\V>=\\V> <W\E = -<W\ (2) 

p a 

and Zjy is a normalization constant. (For the dynamical matrix ansatz see [21].) It 
can be shown that the representations of the algebra (2) are either one dimensional or 
infinite dimensional. The one dimensional representations are the ones which give the 
steady state on the ring( where the matrix element (1) is replaced by a trace), and on the 
open chain when a + f3 = 1. The infinite dimensional representations are used for the 
open chain when a + ft 7^ 1. Although the details of the one dimensional representations 
are important for the properties of the steady state in the open chain, for the periodic 
chain where the number of particles are fixed, once the existence of a one dimensional 
representation is shown, it follows that all configurations are of equal weights and cal- 
culation of all equal time quantities can be done simply by using combinatorics. In [22] 
it has been shown that the MPA is really not an ansatz and in fact the steady state 
of any one-dimensional stochastic process governed by a Hamiltonian with only nearest 
neighbor interaction can be put in the form (1). In the general case where each variable 
r k takes p + 1 values say 0, 1, ...p, one writes again the steady state, in the same form as 
in (1). The p+1 operators Di := -D(rj) and the vectors \V > and < W\ should satisfy a 
set of algebraic relations [22] written compactly as follows: 

h B {V®V) = X®V-V®X (3) 

(h N V- X)\V >=0 <W\(h 1 V + X) = (4) 

3 



where h B is the part of the bulk Hamiltonian acting on two neighboring sites, h 1 and 
h N are boundary terms, and (T>)i = Di and X are two column matrices. Usually X is a 
suitably chosen numerical matrix. However, finding the representations of the obtained 
algebra is by no means an easy problem and in fact may well be harder than the original 
problem. In this sence the work of Kreb and Sandow [22] implies that the MPA does not 
facilitates much the search for the steady state of one-dimensional stochastic systems or 
the ground state energies of generalized spin chains. However one can reverse the problem 
and turns the observation of [22] into a virtue. Thats, one can postulate a consistent 
algebra with nontrivial representations and then find via (3-4), if a physically meaningful 
process or spin Hamiltonian corresponds to this algebra. In this letter we will go through 
this procedure and define the p-species ASEP. We start from the following algebra 

D t E = -Di + E i = l...p (5) 

Vi 

DjDi = J - J - — j > % (6) 

Vi - Vj 

where the parameteres V{ are finite non-zero real numbers which we order as v± < v-i < 

.... < v p . It is easy to check that this algebra is associative, thats for any k > j > i, 

Dj(DiE) = (DjDi)E and D k (DjDi) = (D k Dj)Di. For details of the derivation and also 

for the representations see [23]. It is now straightforward to check that the following 

Hamiltonian when inserted in (3) gives the above algebra, 

v v 
h B = - v i( E oi ® E i0 - E u <g> E 00 ) - J2(vj - Vi)(Eij ® Eji - E n <g> E u ) (7) 

i=l j>i 

where the choice made for X is as follows: x = —l,Xi = Vi/p, and V = E, (V) i = 
-^-.The factors - are for later convenience. Here the matrices En have the standard 
definition (E^ki = S^kSjj. The local Hilbert space of each site is spanned by the vectors 
|0>,|1 >, ...\p >. The state |0 > means that the corresponding site is vacant and the 
state | j > means that it is occupied by a particle of type j. The Hamiltonian (7) describes 
a process in which particles of type j hop with rate Vj and when they encounter particles of 
type i, with Vi < Vj they interchange their sites, as if fast particles stochastically overtake 



slow particles with a rate Vj — V{. This model seems to be very natural as a simple model 
of one-way traffic flow. Real models of traffic flow are too complex to be amenable to 
analytical treatment. One species ASEP is a very simple model of traffic which has 
been extensively studied in the past few years by analytical methods. Compared with 
f-ASEP, our model,while still based on a nice algebraic structure and hence amenable 
to analytical treatment, incorporates new more realistic features of traffic flow, namely 
the existence of different intrinsic speeds for the cars and the possibility of overtaking 
or passing [24]. Moreover, while 1-ASEP is a suitable model for one-lane traffic flow the 
present model, at least in those situations where the probability of cars riding side by 
side is small, seems to be a good model for multi-lane traffic flow. Note that when fast 
particles with speed of say v 2 reach slow ones with speed v\ < v 2 they pass them only 
stochastically. That is, in a time interval dt a fraction {y 2 — Vi)dt of the fast cars overtake 
and the rest are stopped behind the slow ones. This means that although a fast car 
has an intrinsic speed v 2 > v\ in an empty road,its effective speed depends on both the 
distance and the relative speed of the car ahead. Considering a two particle system ( a 
fast one with coordinate x 2 behind a slow one with coordinat Xi) we can write : 

V 2 Cff -=J t <X2 >= V2 ~ V 1 P(~1, t) (8) 

where P(—l,t) is the probability of the fast particle being one lattice site behind the 
slow particle. The above equation can be obtained after lenghty calculations starting 
from the master equation of the two particle system in coordinate space, although it is 
obvious from the very definiton of the process. It is intuitively clear that P(—l,t) is an 
increasing function of v 2 — v\ and a decreasing function of the initial seperation x\ — x 2 . 
which implies a realistic behaviour for v 2 ^. 

Moreover we have shown elsewhere [23] that the same algebra (5,6) describes p-ASEP on 
an open system where particles of speed V{ enter the left with rate — and leave the right 
with rate Vi + (3 — f . In the open system the unit of time is set so that the average speed 
of particles is unity. Hence a and (3 are the total arrival and average departure rates 
respectively. The forms of these rates and their dependence on speed again is a natural 
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property of highway traffic flow. So much for the relevance of this model to traffic flow. 



2 Correlation Functions 

In the following we restrict ourselves to the periodic boundary condition and calculate 
the steady state and some equal-time correlation functions. 

The algebra (5,6) has a one parameter family of one dimensional representations, namely 
E — -, Di = - J5 — where a is a free parameter. This means that in the steady state 
all configurations have equal weights. If there are rrii particles of type i, and the total 
number of particles is M = mi + rri2 + ...m p ,then the probability of all the configurations 
are equal to : 

C = -^m 1 \m 2 \...rn p \(N - M)\ (9) 
We define the following functions: 



n^{r k ) 
n{r k ) 
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= n ^(r k )+n^(r k ) + ...n^(r k ) 
= v 1 n^(r k )+v 2 n^(r k ) + ...v p n M (r k ) (fO) 



The one point functions < n~ k >, < n k > and < v k > of these quantities give respectively 
the average number density of particles of type (i),the average number density of all types 
of particles and the average intrinsic speed of particles at site k. Due to the uniformity 
of the measure none of the correlation functions will depend on the site indices. It is 
simple to show the following: 

!m,(mj — 1) • ■ 

Higher correlation functions have similar forms (i.e < n^n^nf* >= ^^Zlf^J^ )• ^ s 
a sample calculation we note that: 

<nf>=<nf>=Y,&-n>P(T 1 ,...,T N )=C ]T 1 (12) 

{r a } V2,T3, — T N 



The last sum is the number of ways the rest of particles (except (i)) can be distributed 
on the ring. Combining this with the value of C given in (9) gives the result. From (11) 
the following quantities are calculated: 

M M M(M-l) , , 

<n q >=- <v q >=-<v> <n q n r >= N{N _ 1) (13) 

M 2 < v > 2 -M < v 2 > 

< V- >= WiWrrT) < 14 > 

Here the averages < v >,and < v 2 > are taken with respect to the population present 
in the system, i.e: 

m\V\ + m 2 V2 + ...m p v p 

< v > ' = 77 

M 

while the left hand side averages are taken with respect to the steady state configurations 
of the particles on the ring. The result (14) is the new quantity which can be defined 
in this process and is absent in 1-ASEP. From (13) and (14) one obtains the correlation 
function for intrinsic speeds ( or in fact types of particles ) at different sites: 

g(p) :=< v k vi > - < v k >< vi >= JT^jip < v > 2 - < v 2 >) (15) 



where p = ^ is the density . It is seen that the maximum correlation exists at p = and 
P = an d the minimum correlation at p = 2<l>z AH the other correlation functions 
can be calculated in this manner. For example some lenghty calculation will give 

M 3 < v > 3 -3M 2 < v >< v 2 > +2M < v 3 > 
< W >= N(N-l)(N-2) (16) 

3 Boundary Induced Negative Currents 

An intersting physical quantity to consider is the average current of each type of particles. 
From the form of the process one can write the following continuity equation for the 
density of each type of particles. 

± < n ii) >_< jii) > _ < jii) > 

^ ~ n k > — ^ J k > ^ J k+1 > 



where the current of particles of type (i) is given by: 

< jW >= vi < 4le k > + 5> - «;) < nf_ in f > - J2(v 3 - v t ) < «W > (17) 

j<i j>i 

and e k = 1 — X^n^. 111 fact < n k-i e k > is the probability of site k-1 being filled with an 
(i) particle and site k being empty. . From (14,17) the average current is calculated to 
be: 

<J^>= w ^- Y) {Nv % -M<v>) (18) 
Finally we obtain the total current as : 

< J >: = < J« + J<» + ...J*> >= < » > (19) 

again in accord with the 1-ASEP result. 

The interesting point is that for those particles whose hopping rates are less than 
^ < v >, the currents < jW > turns out to be negative. One may think that this result 
is expectable, once we remind of the exchange of particles which takes place in the process 
i.e; particles ordinarily hop forward but when encountered by faster particles from the 
left they hop backward. However negative currents are induced only on the ring and it 
has been proved rigourously [23] that in an open system all the currents are proportioanl 
to one single current. This is an example of how in a nonequilibrium situation boundary 
conditions drastically affect the behaviour in the bulk. A simple way for controlling 
negative currents by varying the density is as follows: As far as M < M c := N ^™ n where 
v min is the lowest speed of the particles there is no negative current in the system. One 
can now add particles of speed < v > to the system which increases only M without 
changing < v > . For all M > M c negative currents will be developed in the system. 
We would like to stress that the merits of this inherent partial asymmetry are quite 
remarkable. On the one hand, one is tempted to think that this model may also have 
some relevance as a two-way traffic model and on the other hand, the algebra that 
describes this process is much simpler than the commutator-like algebras which one may 
try to define for a multi-species partialy asymmetric process from the begining. In fact 
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the attempted algebras of this later type [25,26] do not have even a general and elegant 
form for arbitrary p. 



4 A note on Shock Waves 

For the following discussion we consider open boundary conditions. In this case the one 
dimensional representations give the steady state for the case a + ft — 1 [23] . The steady 
state is a translationaly invariant Bernouli measure in which the density of particles is 
constant along the chain. In the p — > oo limit, where the speeds of the particles are 
taken from a distribution P(v), the total density and the total currents are given by [23]: 

M - x|& '<»> - xtW (20) 

where 

A(a) = f° —H—P( v )dv (21) 
J a v — a 

The parameter a > when eliminated between the above expressions gives J as a function 
of p. We will show that for a large class of distributions, J(p) is a convex function. The 

only condition that P(v) should satisfy is that lim„ > a ^"^ 2 = 0. To prove this we 

rewrite (20) as: 

f°° v 

J — ail - p) . - . = a / P(v)dv (22) 

(1 — p) Ja v — a 

from which = 0(1 - p) - 2f p . Denoting the integral / Q °° ^P(v)dv by I k ; noting 
that -jj^Ik — klk+i for k — 1, 2 and differentiating the second relation of (22) twice with 
respect to p, we find after some algebra 

dp 2 (l-p) 3 (/i + «/ 2 ) 2 1 ) 

which is always negative. This then hints and only hints at the possiblity of shock wave 
solutions for the equaion ff + = ff + J'(p)fa = 0> ^ ne simplest of which are in 
the form of step functions with density p_ at the left and p + at the right of the front 
which itself moves with the speed V s h oc k = 3 p~-p + ■ ^ n conclusion, the approach presented 



in this letter, that is, begining from consistent algebras and then finding the process 
may be further pursued to study other one dimensional reaction diffusion processes. It 
may also be possible to map this model to more complex interface growth models. In 
this letter we have been lucky to postulate an algebra which while being simple and 
beautiful, approximately models traffic, an every-day life phenomenon. 
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